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Recent results on the QCD analysis of bound states of heavy qq quarks are reviewed, paying attention 
to what can be derived from the theory with a reasonable degree of rigour. We report a calculation of 
be bound states; a very precise evaluation of 6, c quark masses from quarkonium spectrum; the NNLO 
evaluation of T e^e~ ; and a discussion of power corrections. For the b quark pole mass we get, including 
0(m^/TO^) and 0{al log as) corrections, mf, = 5.020±0.058 GeV; and for the MS mass the result, correct 
to 0{al), 0{m^/m1), fhb{fhi,) = 4.286 ±0.036 GeV. For the decay T e+e", higher corrections are too 
large to permit a reliable calculation, but we can predict a toponium width of 13 ± 1 keV. 



1. INTRODUCTION 

In the present note we are going to review some 
aspects of the QCD analysis of heavy quarkonia, 
tt, cc, 6c and especially bb states. The validity of 
ab initio QCD calculations varies from one case to 
another. For the energy levels of toponium with 
n < 5, n being the principal quantum number, and 
for the energy of the ground state in bottomium, 
we have a very satisfactory situation, almost com- 
parable to that in positronium calculations. For 
bottomium hyperfine splitting, and for the energy 
of the charmonium ground state, the situation is 
less favourable, although reasonably under con- 
trol. Less favourable still is what one has for the 
energy levels of bottomium with n = 2. Here NLO 
and NNLO perturbative and the leading nonper- 
turbative corrections are comparable to the LO 
results. The same is true for the decay T e"'"e~, 
and even for toponium decay, T — > e+e^, higher 
corrections are a bit too large for comfort. For 
higher states, for the decays of n = 2 bottomium 
or for the ground state of charmonium, a rigorous 
QCD evaluation is out of the question. There are 
methods that have been devised in the literature 
to deal with this, essentially more or less justified 
models; we do not discuss them now. The inter- 
ested reader may find information, and references, 
in e.g. my Lisbon Icctures^^l. 

This does not mean that even in the favourable 
cases all questions are settled, and in particular we 
discuss here power-like corrections, at short dis- 
tances, where there is still some controversy. 



2. NONRELATIVISTIC QCD POTENTIAL, 
AND CORRECTIONS 

To analyze the lowest states of heavy quarks one 
proceeds as follows: first, in the nonrelativistic 
(NR) approximation one finds the potential that 
governs the dynamics. This potential is evalu- 
ated to increasing orders of perturbation theory; 
at present it is known to one and twol^l loops. One 
can thus write 

i/ = i/(o)+i/i„g; (2.1a) 

where H'-°^ may be solved exactly and contains the 
coulombic part of the interaction. For qq states, 

ff(0)=2m+^A-^^M^, (2.1b) 
m r 

a.(,.^)=a.(,.^)|l + cW^+c(^)4|; 



(2) ( « /tT^ ^ 

' =7E ( ai/3o + y) + (Y2^'^jT ^^' 

The expression for the ai may be found in e.g. 
ref. 1, m is the pole mass of the quark, and i/iog 
is 

27r r ^ 4^2 ^ 

, -Cpal ( fD .f^i . lBf3o \ logr^f 



* Talk given at the conference " Hyperons, Charm and Beauty Hadrons", Valencia, June 2000. To be published 
in Nucl. Phys. Suppl. 



One solves exactly the coulombic Schrodinger 
equation, 

Hmr,i=E^°^^>I'^i, (2.2) 

treating the terms in Hiog as perturbations. Rel- 
ativistic corrections arc then addc;(l. to be also 
treated as perturbations. These relativistic correc- 
tions include spin-dependent ones, known at tree 
level since a long time (they are essentially like for 
positronium) and to one loop from refs. 3,4. More- 
over, and because of the nonabelian character of 
QCD, we have mixed relativistic -one loop correc- 
tions for the spin-independent piece, first evalu- 
ated in ref. 5 using the method of equivalent po- 
tentials of Gupta et al. (see e.g. ref. 4). Concen- 
trating on the spin-independent hamiltonian we 
then write 



rrc.p. _ rr(0) 

-"si - 



log" 



-1 

4m3 



2mr^ 
(2.3) 

and bi = \Cf — '2'Ca- Here the superindex "e.p." 
means that the hamiltonian was obtained with the 
method of equivalent potentials. 

The hamiltonian has been re-evaluated recently 
by Brambilla et al.^^l; these authors generalize 
the nonrelativistic (heavy quark) effective theory 
(NRQCD^^l), matching then the Green's functions 
to a potential description (pNRQCDl^l). In this 
way they find. 



pHQET 
SI 



-1 



4m'^ 



H —^{r) + - . 



2r3 



Cfo, 1 

2to2 '■r' (2.4) 



2mr2 

This is different from (2.3). The difference is 
due to the fact that in the derivation of (2.3) one 
works with the S-matrix, while for (2.4) one uses 
Green's functions. The difference between (2.3) 
and (2.4) vanishes when taking expectation val- 
ues between coulombic wave functions, say, solu- 
tions of (2.2): so they will produce the same en- 
ergy spectrum, at least to order af. At higher 
orders, (2.3) and (2.4) would produce different re- 
sults; but also new terms will appear in the two 
formalisms: their equivalence, or non-equivalence, 
has not been proved. For numerical results to 
0{af) for all states see ref. 9; here we present those 



for the ground state energy. We write 



E^Q- = 2m - 



m 



— j \-dp.tEi 



(2.5a) 



The label "p.t." in i?^;*' indicates that we have as 
yet only used results deduced from perturbation 
theory; the full expression would be 

Enl = -E'^j*' + STSipEnl. 

6-!sipEni embodies the nonperturbative contribu- 
tions, to be discussed below. The Sp,t.Eio is, with 

a = 2/mCFOis, 



^ F - ^0 



+ 



[log- 



+ 



3 + 37I - TT^ + 6C3 

12 



(2.5b) 



= ai/3o + i/3i + i7E/3o- 



The (nominally) leading, log ag corrections of 
next order are also known. They include a loga- 
rithmic correction to the static potentially*^!, and a 
relativistic one-loop correctionl^'^^I. With the full 
result (as given in ref. 6) one has, for /x = 2/a, 

S[a^-ioga,]Eio = -m[CF + | Ca] C|. Qf^ (log a;s)/7r. 

(2.5c) 

In the calculations we will include, for the b quark 
case, a correctionl^^l of order ml /ml, 



^lm^/ml]Eio = 2mb 



STfUs ml 
27r m? " 



(2.5d) 



We can invert (2.5) to obtain a precise determina- 
tion of c, b quark masses from those of the J/ip, T 
particles. The determination includes nonpertur- 
bative effects, to be discussed below. As input 
parameters we take the recent determination'^^l , 

A{nf = 4, three loops) = 0.283 ± 0.035 GeV 
[ a,(M|) ~ 0.117 ±0.024] , 

and for the gluon condensate the value (asG^) = 
0.06 ± 0.02 GeV". 

Then comes the matter of the renormalization 
point, /i. As Eqs. (2.5b,c) show, logarithms are 
avoided until order if choosing fi = hq = 2/a, 



and this will bo our central value here. Then, for 

the b quark, we find 

mf, =5.020 ± 0.043 (yl) T 0.005 ((a^G^)) 
(vary ^2 by 25%) 
±0.006 ± 0.017 (other th. uncert.) 
=5.020 ± 0.058 GeV . 

(2.6a) 

To obtain this, we have included the correction of 
Eq. (2.5c), and for the "theoretical error" also that 
estimated in ref. 13 for the pcrturbative evalua- 
tions. The value of the MS mass that corresponds 
to this is, taking into account 0{al) and leading 
0{m'^/ml) corrections!^'*!. 



mbifhb) =4272 ±43 MeV. 



(2.6b) 



We present a summary of QCD calculations of the 
quark masses, with increasing accuracy. 



Ref. 


m6(pole) 


mb{ml) 


mc(pole) 


TY 
PY 
Here 


4971 ± 72 
5065 ± 60 
5022 ± 58 


4401^35 
4455t^^ 
4272 ± 43 


1585 ±20* 
1866i?i^ 



b, c quark masses. (*) Systematic errors not included. 

TY: Titard and Yndurain'^^l . O(a^) plus 0(aj)v, 
0{v^) for m; 0(0?) for fh. Rescaled for A{nf = 
4) = 283 MeV. 

PY: Pineda and Yndurainl^l. Pull 0{ai) for m; 
0{al) for m. Rescaled for A{nf = 4) = 283 MeV. 
Here: This calculation. 0(0,), OiasnvHiTib) ^^'^ 
0{a^s logQs) for m; 0{a^s) and 0{ai''m'i/m1) for rh. 
Values not given for the MS c quark mass, as the 
higher order terms are as large as the leading ones. 

It is to be noted that pure QCD quarkonium 
determinations of m;, have remained remarkably 
stable since the earlier evaluations (ref. 5). [The 
variation of the published values is due to the vari- 
ation of the favoured central value of A, from 200 
to 283 MeV]. Also, the results depend very lit- 
tle on the value of /i chosen. Thus, even taking 
/X = m ~ 5 GeV one gets 

nib = 4901 ± 60 MeV, 

reasonably close to (2.6a), where we had elected 
/i = 2/a ~ 2.88 GeV. The estimates of fhb{fhb) 
are less stable, doubtlessly because of the large 
size of the O(a^) corrections; but it is interesting 



that this more precise spectroscopic determination 
of MS mass of the b quark mass given here agrees 
with recent determinations, based on sum rules, 
also accurate to 0{a^). These give'^^l 

ffib{rfib) = 4260 ± 100 MeV . 

We will discuss MS masses further later on. 

The hypcrfinc splitting can likewise be evalu- 
ated, getting a prediction for the r^b mass: 

M{T) - M{r]b) = 53.3 ± 12.5 MeV . 

The methods discussed for quarkonium can be 
extended with some work to hybrid states like be. 
The advantage here is that, since the masses of 6, c 
quarks can be expressed in terms of those of the 
T, J I tp particles, the theoretical errors are smaller 
than what one would get for charmonium. A re- 
cent 0{ag) calculation!^^! gives 

Mw{bc) = 6323 ± 10 ± 20 MeV, 

the first error being pcrturbative, the second com- 
ing from estimated nonperturbative effects. 

3. POWER CORRECTIONS. THE NONPER- 
TURBATIVE VACUUM; RENORMALONS; 
SATURATION 

As stated above, a calculation such as that in (2.6) 
includes leading nonperturbative effects. These 
are obtained by realizing that the motion of the 
qq pair takes place in the physical vacuum, chock 
full of soft gluons and light quark pairs so that, in 
particular, we have a nonzero value for the gluon 
condensate: 

(vac| : a,G^{0) : |vac) = {a.C^) ^ 0. 

The effects of this were first considered by Leut- 

wyler and Voloshinl^^l (sec also refs.. 18, 5) and 
amount to a shift for (say) the ground state en- 
ergy of quarkonium of 

ne^oia. : ^ (3.1) 



•'NP 



{mCFOis 



eio — 1.5. This is zero to all orders of perturbation 
theory, because so is the gluon condensate. On 
dimensional grounds one expects {agC^) ~ yl^, so 
that (3.1) is of order A^/m^a*. 

This is not the only power correction that may 
appear. We may have renormalons. Let us con- 
sider a one-gluon exchange diagram, for qq scatter- 
ing. If we dress the gluon propagator with loops 



then the corresponding potential, in momentum 
space, is 



V{k) 



477 



fc2 /Jolog(fcVyl2)' 



(3.2) 



and we have substituted the one-loop expression 
for asik"^)- (3.2) is undefined for soft gluons, with 
f,2 ^ follows from the general theory of 

singular functions, the ambiguity is of the form 
C(5(fc^ — yl^): upon Fourier transformation this pro- 
duces an ambiguity in the x-space potential of 
dV{r) = c[sinylr]/r. At short distances we ex- 
pand this in powers of r and find'^^'^^l 



.V{r) ~ Co + Cir^ + . . . . 



(3.3) 



The same result may be obtained with the more 
traditional method of Borel transforms. This co- 
incides with the short distance behaviour of the 
nonperturbative potential as determined by Dosch 
and Simonoyl^^l. 

The situation just described applies for states 

qq at short distances; but not so short that zero 
frequency gluons cannot separate the pair. If this 
last is the case, soft gluons do not resolve the qq 
pair and only see a dipole. The generated renor- 
malon may then be seen (ref. 20) to correspond 
to the contribution of the gluon condensate in the 
Leutwyler-Voloshin mechanism. 

What happens to the Aglietti and Ligeti rcnor- 
malon? In fact, we have other long-distance power 
corrections. It is clear that the pole mass is defined 
purely in perturbation theory, and indeed one can 
check that a renormalon ambiguity appears al- 
ready at one loopP^l. Likewise, the coulombic 
potential is defined so that it vanishes at infinity: 
but, for confined quarks, "infinity" is equivalent to 
the confinement radius, ~ 1/A. Actually, mass 
and potential renormalons cancel for the constant 
and quadratic terms in (3.2), in the very short dis- 
tance regime; this has been verified by a detailed 
calculation in ref. 22, thus substantiating the in- 
tuitive arguments of ref. 20, already mentioned. 



X, 



FIGURE 1. "Filled in" gluon loop. 

Renormalons and the mechanism of Leutwylcr - 
Voloshin are not the only possible sources of power 
corrections; we next consider saturation. We note 
that the ambiguities we have found arc associ- 
ated with small momenta or, equivalently, long 
distances. However, the singularities are clearly 
spurious. Indeed, not only the theory should be 
well defined but, because of confinement, long dis- 
tances are never attained: the theory possesses an 
internal infrared cut-off of the order of the con- 
finement radius, R ~ A~^. To try and implement 
it we consider the gluon propagator. To one loop 
it gets a correction involving the vacuum polariza- 
tion tensor. Neglecting quarks this is, in x-space, 
given by an expression like 

+ ■■■ . 

We can take into account the long distance inter- 
actions by introducing a string between the field 
products at finite distances. In matrix notation 
for the gluonic fields, B'' = taB^, this is incorpo- 
rated by replacing 

^ S«(yi)P (^exp ij^' dz^^B^{z)^ B^iy^). 

The process may be described as "filling the 
loop" (see Fig. 1) by introducing all exchanges be- 
tween the gluonic lines there. If we furthermore 
replace the perturbative vacuum |0) by the non- 
perturbative one I vac), then we get a dressed prop- 
agator 



An 



/3olog(M2-Kfc2)/yl2' 



and is related to the gluon condensate at finite 



distances, (G(.T)G'(0))vac- (for more details and 
references, see e.g. Simonov's Iccturcs'^^l). 

This suggests a saturation property of the cou- 
phng constant at small momenta (long distances) 
so that the expression for the running coupling 
constant should be modified according to 



/3ologfcVyl2 
.af (fc^) 



(3.4) 



/?0log(P +M2)/vl2' 



It is certain that an expression such as (3.4) in- 
corporates, to some extent, long distance proper- 
tics of the QCD interaction. For example, if we 
take (3.4) with M = yl in the tree level poten- 
tial for heavy quarks, this becomes the Richardson 
potentiaP^^ 



FW(k) 



_ 47rCFa,(k2) 



/Jok2log(fc2+yl2)/yl2 



When one has k2 ^ A^, the short distance coulom- 
bic potential is, of course, recovered. For k2 <C A^, 
however, 

y(0),sat.k) ^ I&tt^C'fA^ 

whose Fourier transform gives 

y(0),sat^^^ _ (constant) x r, 

i.e., a linear potential. Indeed, a reasonably accu- 
rate description of spin-independent splittings in 
quarkonia states is obtained with such a poten- 
tial. However, this long-distance linear potential 
induced by saturation in the Richardson model is 
the fourth component of a Lorentz vector, while 
we know that the Wilson linear potential, as ob- 
tained, e. g., in the stochastic vacuum model or 
in lattice calculations, should be a Lorentz scalar. 
This is also obvious on phenomenological grounds 
as this potential has to provide attraction for 
quaric-quarJc states in e.g. baryons. It thus follows 
that the linear potential obtained from saturation 
can be only of limited phenomenological use. 

Saturation also gives a linear correction to the 
short distance coulombic interaction. The possi- 
bility of such a correction has been discussed by 



several people; sec, for example rcfs.. 20, 23. Per- 
sonally I am not impressed by these arguments. 
It is clear that the QCD perturbative series is not 
convergent; in the case of quarkonia, the coeffi- 
cients are not even analytic as one has terms in 
log as (cf. (2.5c) and Sect. 4 below). The results 
one finds are then dependent on how one sums 
and, moreover, it is not guaranteed that results 
valid for large orders of perturbation theory will 
hold in the real world, where one knows at most 
three nontrivial terms. It is the author's belief 
that only if the summation method is rooted on 
solid physics it is likely to represent an improve- 
ment; otherwise, estimates of nonperturbative ef- 
fects become pure guesswork. In this respect, the 
method of taking into account the nonperturba- 
tive nature of the physical vacuum by consider- 
ing the effect of nonzero values for the correlators 
stands some chance of being meaningful, as indeed 
phenomenological calculations seem to indicate. 

4. FURTHER DISCUSSION OF THE MS MASS: 
A MATTER OF CONVERGENCE 

Write, for a heavy quark, 



m(m) =m/{l + 6i + 62 + S3 + ■ ■ ■}, 



(4.1a) 



with the coefficients C„ given in refs. 14 (see also 
ref. 11). Numerically one has 



5i(6) =0.090, (5i(c) =0.137, 
62{b) =0.045, 52(c) =0.108, 
53(b) =0.029; 53(c) =0.125. 



(4.2) 



Thus, the series relating pole and MS masses is 
at the edge of the region of convergence for the 
b quark, and clearly diverges for the third order 
contribution for the c quark. 

We can use (4.2) and the formula for M(T) 
in terms of the pole mass to express the former 
in terms of Wf,: One finds, for the perturbative 
contributions, neglecting m^/ml, and with A as 
before, 

Q!s(m) 



M(r) =2mb{mb)\l + CF^ 



+7.559 



as (to) 



-I- 43.502 



(4.3) 

(We have taken nj = 4). This does not look 
particularly convergent, and is certainly not an 



improvement over the perturbative expression us- 
ing the pole mass where one has, for the choice 
ttsin = Cprribas) and still neglecting the masses 
of quarks lighter than the 6, 

M{r) =2mf,{l - 2.193 (^)^ - 24.725 (^)^ 

-458.28 (— )^ + 897.93 [loga,] (— }• 

(4.4) 

To order three, (4.4) is actually better than (4.3); 
thus, and at least to third order, it is unclear that 
there is a connection between the rate of conver- 
gence and the use of pole or MS mass. What is 
more, at least the size of the known power cor- 
rections (namely, those stemming from the gluon 
condensate) do not favour the expression in terms 
of the MS mass. Indeed, if we evaluate the gluon 
condensate corrections to the direct formula (4.3) 
for M{T) in terms of fhbijhb), it turns out that 
they are larger than than what one would have 
for the expression in terms of the pole mass, (4.4) 
(~ 80 against ~ 9 MeV), because of the definition 
of the renormalization point. 

5. DECAYS OF QUARKONIA 

Besides energies of bound states, once one has an 
effective interaction it is possible to evaluate, at 
least in principle, decay rates such as T — » e^e~ . 
To relative order as (NLO) this has been done 
in ref. 5 (see also ref. 9) using the method of 
variations and the evaluation of ref. 25 for the 
hard part, ^rad- We give here the exact result of 
a Rayleigh-Schrodinger perturbative calculation. 
One finds. 



r{T ^e+e-) = [1 + <5wf + + <5np]' 



QbCKQED 



M{V) 



2 



c , r I 3/3o /, aM TT 



+- 



3c(i) 




(5.1) 

given in (2.1) and the nonperturbative contri- 



.(1) 



bution is 



3np 



1 [ 270 459 

2 L 108 800 



1838 
2 890 000 J 



7811 71"(asG2) 



The corrections are very large. Because of this 
the calculation is unreliable, even for bb, and fails 
completely for cc. One might hope that this would 
be arranged by the 0{al) corrections (NNLO); 
but this appears not to be the case. We have, first, 
"hard" corrections'^^l, similar to 5rad; and correc- 
tions to the wave function, given in refs. 11,27 
where we send for the (rather lengthy) explicit for- 
mulas. 

Numerically, and with the width in keV, 





e+e-) 


LO 


NLO 


NNLO 




: m 


0.41 


1.22 


1.13 




2/a 


0.73 


0.55 


0.80 



The experimental figure is 

r{T ^ e+e-) = 1.32 ± 0.04 keV . 

Clearly, the large NLO and NNLO perturbative 
corrections, both of similar size, and of the leading 
NP correction, make the theoretical result unsta- 
ble, as the Table shows clearly. 

For the (perhaps measurable) toponium (T) 
width we get, for rrit = 175 GeV and still in keV, 



r{T e+e-) 


LO 


NLO 


NNLO 


p = mt 
fi = 2/a 


6.86 
10.24 


10.53 
10.91 


13.0 
13.5 



The situation has improved with respect to 
what we had for bottomium, but there is still a no- 
ticeable dependence on the renormalization point, 
and on the order of perturbation theory consid- 
ered. We conclude on an estimate of some 11 — 14 
keV for the width. 
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